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BREAKING  WAVE  SPECTRUM  IN  WATER  OF  FINITE  DEPTH  IN  THE 


PRESENCE  OF  CURRENT 


PART  I:  INTRODUCTION 

1.  There  are  many  forms  of  wave  energy  spectrum.  All  of  these  spectra, 
however,  are  for  specific  conditions.  For  example,  the  Pierson-Moskowitz 
spectrum  is  for  a  fully  developed  sea,  the  Joint  North  Sea  Wave  Project 
spectrum  is  for  a  f e tch- 1 imited  developing  sea,  and  the  Wallops  spectrum 
(Huang  et  al.  1981)  is  derived  based  on  wave  dynamics  but  without  considering 
wave  breaking. 

2.  When  conditions  differ  from  those  for  which  these  spectra  arc  in¬ 
tended  or,  as  the  waves  move  into  regions  where  the  conditions  are  changed, 
these  spectra  undergo  corresponding  changes.  For  example,  as  the  steepness  of 
the  wave  increases,  wave  breaking  occurs;  the  Wallops  spectrum,  which  does  not 
consider  wave  breaking,  must  be  modified.  This  problem  was  treated  recently 
by  Yuan,  Tung,  and  Huang  ("19861  and  by  Tung  and  Huang  (1987)  for  deepwater 
waves . 

3.  As  the  waves  propagate  from  deep  to  shallow  water,  wave  breaking 
takes  place  when  they  reach  the  surf  zone.  There  have  been  a  number  of  publi¬ 
cations  on  the  subject  of  wave  breaking  in  shallow  water  such  as  those  by 
Battjes  and  Janssen  (1978)  and  Thornton  and  Guza  (1983).  These  authors  used 
the  energy  flux  balance  equation  including  energy  dissipation,  and  the  results 
are  in  good  agreement  with  measurements.  The  equation,  however,  must  be  inte¬ 
grated  numerically,  and  the  methods  do  not  give  the  breaking  wave  spectrum 
direct  1 y . 

4.  Similarly,  when  a  wave  train  encounters  ar  adverse  current,  wave 
breaking  results.  The  method  usually  employed  to  obtain  the  spectrum  of  the 
waves  interacting  with  current  is  to  first  resort  to  the  classical  energy  flux 
balance  equation  without  considering  wave  breaking  (Huang  et  al.  1972  and 
Hedges,  Burrows,  and  Mason  1979).  To  account  for  the  effect  of  wave  breaking 
‘■n  the  wave  spectrum.  Hedges,  Burrows,  and  Mason  tl979)  applied  the  equilib¬ 
rium  range  spectrum  to  limit  the  spectral  ordinates.  The  equilibrium  range 
spectrum,  however,  contains  a  numerical  constant  whose  value  is  difficult  to 
specify.  Furthermore,  the  equilibrium  range  spectrum  on  1 v  applies  to 


! 


frequencies  much  higher  than  those  corresponding  to  the  peak  of  the  spectrum 
and  therefore  cannot  be  extended  to  cover  the  range  of  frequencies  where  most 
of  the  wave  energy  resides. 

5.  In  this  study,  we  extend  the  method  introduced  earlier  (Yuan,  Tung, 
and  Huang  1986,  Tung  and  Huang  1987)  for  a  deepwater  breaking  wave  spectrum  to 
waves  in  water  of  finite  depth  and  apply  the  results  to  the  situation  where 
the  waves  encounter  a  current.  The  method  consists  essentially  of  first 
assuming  that  there  exists  an  original  ideal  wave  train  at  the  locale  under 
consideration,  the  spectrum  of  which  is  obtained  from  the  equation  of  energy 
flux  balance  without  considering  wave  breaking.  Py  imposing  the  Miche  wave 
hreaking  criterion  (Battjes  1  974)  ,  an  expression  for  the  elevation  of  the 
breaking  waves  is  established  in  terms  of  the  original  ideal  wave  elevation 
and  its  second  derivative  which  are  assumed  to  be  jointly  Gaussian.  Based  on 
this  breaking  wave  model,  the  expressions  for  the  mean  value,  the  mean-square 
value,  and  the  spectrum  of  the  breaking  waves  are  derived.  These  results  arc 
then  applied  to  the  case  in  which  a  unidirectional  deepwater  wave  train, 
propagating  normally  toward  a  straight  shoreline  over  a  gently  varying  sea 
bottom  with  straight  and  parallel  contours,  meets  a  steady  current  whose  flow 
velocity  is  uniformly  distributed  in  the  vertical  direction.  Numerical 
results  are  obtained  and  given  in  graphical  form.  The  simpler  breaking  wave 
model  for  deepwater  waves  is  first  presented  and  modified  for  waves  in  water 
of  finite  depth. 

6.  It  is  emphasized  here  that  the  studies  carried  out  in  this  report 
are  based  on  heuristic  wave  breaking  models  and.  simplified  current  and  coast 
conf i gurat ions.  A  number  of  approximations  are  introduced  in  the  derivations, 
but  the  results  have  not  yet  been  checked  against  either  field  or  laboratory 
experiments.  It  is  clear  that  the  models  have  yet  to  he  modified  and  that 
more  detailed  studies  should  be  performed  to  examine  the  effect  of  utilizing 
various  spectral  forms  for  the  original  ideal  waves. 
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PART  II:  BREAKING  WAVE  MODELS 


7.  Stokes  (1880)  showed  that  in  deep  water,  when  the  vertical  downward 
acceleration  at  the  crest  of  the  wave  reaches  a  value  of  0.5g*  (g  being  gravi¬ 
tational  acceleration),  the  wave  breaks  and  its  amplitude  is  reduced  according 
to  the  ratio  of  0.5g  and  the  magnitude  of  the  acceleration  of  the  original 
ideal  wave  at  the  crest.  The  following  equation  expresses  this  relationship: 


a 


b 


0^  = 

2  2 
aw  w 


n) 


where 

a^  =  amplitude  of  the  breaking  wave 
a  =  amplitude  of  the  ideal  wave 
v  =  frequency  of  this  ideal  wave 

8.  I.onguet-Higgins  (196°)  applied  this  criterion  to  a  narrow-band  wave 
train  in  which  the  amplitude  of  the  breaking  wave  is  given  by 
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where 


/  JS(  .  )c! . 
sTTTdT"- 

is  the  characterist ic  wave  frequence  and  SO.) 
ideal  waves. 


i  n 


is  the  energv  spectrum  of  the 


v.  T('  obtain  the  spectrum  of  the  breaking  waves,  we  assume  'Phillips 
l'W'  that  the  wave  breaks  whenever  the  local  vertical  downward  acceleration 
at  any  point  on  the  surinte  reaches  a  fraction  of  the  gravitational  accelera¬ 
tion.  Keierring  to  Figure  I,  let  T ( t 1  and  '  (t)  represent,  respectively, 

the  elevation--  ,-t  the  ideal  and  breaking  waves  at  a  fixed  point  in  space  where 


F-  - 1  c  rivei.i  nee ,  svmhols  and  abbreviations  are  listed  in  the  Notation 

l  Append  i  >  B  i  . 
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t  is  time.  Wave  breaking  may  cake  place  at  points  such  as  A  and  B  where 
: ( t 1  ■  D  .  (Here,  and  hereafter,  overdot  denotes  differentiation  with  respect 
to  time.''  At  points  such  as  A  where  4  ( t )  >0  and  when  r(t)  <  -Kg  (K  i  s 
shown  to  range  between  0.4  and  0.5  (Ochi  and  Tsai  1983)),  the  breaking  wave 
elevation  is  given  by 


(t) 


■;(t) 


-Kg 
'( t'< 


:  41 


This  expression  is  a  restatement  of  liquation  I  or  Equation  2;  that  is,  when 
the  wave  breaks,  the  local  wave  elevation,  is  reduced  according  to  the  ratio  of 
Kg  and  the  magnitude  of  the  local  acceleration  of  tin.  ideal  wave. 

!0.  At  point  H  where  'it'1  •  0  and  when  ' !  t  i  *'  -Kg  ,  the  breaking 
wave  elevation  is  given  by 


( t 


’ft  1 


-fit) 

Kg 


11.  Based  on  the  above  considerations  and  noting  that  no  wave  breaking 
takes  place  when  f(t)  ■  -Kg  ,  in  which  case  :,( t  '  remains  unchanged ,  ' 1  t 

nav  he  written  as 
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■h 


H(-r 


Kg)  H  ( rj 


Kg)  H  (  -  ’ )  +  fH  (  +  Kg' 
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where  H  ( • )  is  the  Heaviside  unit  step  function  and,  for  brevity,  the  argu¬ 
ment  t  in  t,  (t)  ,  c(t)  ,  and  ’ ( t )  is  omitted.  In  Equation  6,  the  first 

b 

and  the  second  terms  correspond  to  the  points  such  as  A  and  B  in  Figure  1  when 
wave  breaking  occurs,  and  the  third  term  merely  states  that  4  remains  un¬ 
changed  as  long  as  s(t)  >  -Kg  regardless  of  the  point  under  consideration. 

12.  In  Equation  6,  the  breaking  wave  elevation  Q.  is  a  nonlinear 

b 

function  of  t  and  4  ,  the  elevation  and  its  second  derivative  of  the 
original  ideal  waves  which  are  assumed  to  he  stationary  and  jointly  Gaussian 
with  zero  mean  values.  The  determination  of  the  mean  value,  mean-square 
value,  and  the  spectrum  of  may  therefore  be  achieved  in  a  straightforward 

manner  f Papou 1  is  1965). 

13.  In  water  of  finite  depth,  for  a  single  wave,  the  breaking  wave 
amplitude  is  (Battjes  1974' 


=  0 . 44d 
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tanh  kd 
kd 


(7) 


where  d  is  the  local  water  depth  and  k  is  the  wave  number.  The  above  may¬ 
be  expressed  ap 
that  is,  using 


be  expressed  approximately  in  terms  of  k^  ,  the  wave  number  in  deep  water; 


k  =  k  tanh  kd 
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and 


k  =  k  Vtanh  k  d 

o  ’  o 

we  h  a  v  e 

"anh  k  d 

a,  =  0.44  - - - — 

b  k 
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(9) 


(10) 


The  deepwater  wave  number  is 


where  w  is  Che  wave  frequency  which  remains  constant  independent  • 
depth  and  is  related  to  the  local  wave  elevation  1  and  the  sur  i  ace 
acceleration  at  the  point  under  consideration  hv 


14.  For  random  waves  in  water  of  finite  depth,  for  reasons  purely  of 

mathematical  convenience,  we  replace  in  Equation  10  by  and  k  it 

the  denominator  bv  k  =  -’/eg  but  substitute  the  same  in  the  numerator  hv 

o 

k  =  ..  “  /  g  ,  the  characteristic  wave  number  in  deep  water,  so  that,  i  or  point 
o 

such  as  A  in  Figure  !,  the  breaking  wave  elevation  t  takes  the  same  torn  "■ 
that  shown  in  Equation  4  where 

K  =  0.44  tanh  k  d  ‘  ‘ 

o 

which  reduces  to  K  =  0.44  in  deep  water. 

]>.  Following  the  same  argument  leading  to  Equation  6  for  deepwater 
waves ,  we  see  that  the  breaking  wave  elevation  in  water  of  finite  depth  is 
also  given  by  Equation  n  with  K.  replaced  by  Equation  I ’* .  In  arriving  at 
this  expression  we  have  ignored  the  situation,  where  the  magnitude  >'•:  the 
negative  wave  elevation  mar.  exceed  the  waiter  depth.  he  chose  not  to  consider 
such  pass  ib  i  1  i  tv  :.d  restrict  the  application  of  tin  node  I  to  regions  that  ar 
not  uraiu I v  shallow  where  other  factors  such  as  bottom  t  r ict ion  man  come  into 
n  1  a.-.  ,  the  line. in  :nd  i-aussi  ar.  ns  sump t  i  ons  of  the  original  waves  are  no  longer 
valid,  and  the  mail.;  may  be  o’.erlv  strained. 

Is.  !  n  subsequent  derivations  of  the  expressions  for  the  mean  value, 
t  he  mean-square  value,  and  the  spectrum  of  ‘  ,  the  second  term  in  Equation 

is  ignored  based  on  the  cons i de ra t i on  that  the  probability  of  occurrence  of 
negative  peaks  such  as  point  B  in  Figure  1  is  usually  small,  especially  when 
the  spectrum  of  the  waves  under  cons i derat i on  is  reasonably  narrow.  In,  this 
wav,  the  derivation  is  much  shortened,  and  our  computation  shows  that  the 
error  incurred  hy  ignoring  the  cert  r.d  term  in  Equation  f>  is  indeed  impercepti 


PART  III:  MF AN  VALUE,  MEAN-SQUARE  VALUE,  AND  SPECTRUM  OK  r 

h 

17.  Although  the  original  ideal  waves  are  a  zero  mean  process,  from 
Equation  6  it  is  obvious  that  the  elevation  of  the  breaking  wave  is  not.  From 
Equation  6  (with  the  second  term  deleted),  it  is  not  difficult  to  show  that 
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where 

Ef  •] 
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;:(•) 

I.e  t 


expected  value  of  the  quantity  enclosed  in  brackets 
wave  breaking  parameter  (defined  in  equation  19) 
spectral  bandwidth  parameter  (defined  in  Equation  18) 
probability  function  (defined  in  Equation  20) 
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where  S  (  .. ) 


is  the  spectrum  of  the  original  ideal  waves.  The  quantity 
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(4) 

r  r 


(18) 


lies  between  zero  and  unity  and  is  known  as  the  bandwidth  parameter  of  S(w) 
(Cartwright  and  Longuet-Hi ggins  195b)  and 
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,  ,  KS 

nTT1* 

is  n  measure  oi‘  the  extent  of  wave  breaking  as  will  be  shown  later, 
f  unct ions 
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where 


P 

0 


probability  function 

parameters 

parameter 

probability  function 


0  <  x  i 


j  Z  (vldy 


are  probability  functions  ( Abramowi t z  and  Stegun  19681,  and 


!,fxl 


is  the  exponential  integral  (Abramowi t z  and  Stegun  1 9h8 ) ,  and  v  is 

variable.  It  is  seen  that  the  mean  value  of  r,  is  a  nonzero  const 

n 

depends  on  the  value  of  water  depth  d  through  K  in  .•  in  Kquati 
the  zeroth,  second,  and  fourth  spectral  moments  ot  the  ideal  waves. 


(19) 

The 


(20) 


(21) 


(2  21 


i  1 1 

a  dummv 
ant  which 
on  i  't  and 
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19.  To  obtain  the  spectrum  of  £,  ,  we  first  form  its  autocorrelation 

b 

function.  For  convenience,  let  subscripts  1  and  2  refer  to  quantities 
evaluated  at  time  instants  tj  =  t  +  i  (where  r  is  time  lag)  and  t^  =  t  , 
respectively.  Furthermore,  let  H  =  H(r)  ,  H"  =  H(t;  +  Kg)  ,  and 
H"  =  H(-'  -  Kg)  .  By  anticipating  that  is  stationary,  the  auto¬ 

correlation  function  of  r  ,  denoted  P,  (t)  ,  is,  from  Equation  6  (with  the 

b  b 

second  term  deleted) 


The  expected  values  in  Equation  26  involve  the  random  variables  Cj  >  Co  ’ 
r  ,  and  which  are  jointly  Gaussian  with  zero  mean  values.  These 

expected  values  mav  all  he  obtained,  although  the  task  is  tedious.  In  Ap¬ 
pendix  A,  the  last  expected  value  of  Equation  26  is  evaluated  to  illustrate 
the  techniques  employed  in  obtaining  these  expected  values. 


20.  The  resulting  autocorrelation  function  is  a  nonlinear  function  of 

(2) 

the  correlation  functions  r^Ci)  =  *  rl2  ^  =  *  and 

(4) 

r^  (t)  =  E[^£0]  of  the  original  wave  elevation  t,  and  its  second 
derivative  t.  evaluated  at  time  instants  tj  and  t?  .  The  autocorrelation 
function  ( t )  ,  viewed  as  a  function  of  the  correlation  coefficient 

functions  r^9(x)/r  ,  r^CiO/r^  ,  and  rj^  (t) /r^  may  be  expanded  by 

Taylor's  series  (Borgman  1965).  By  retaining  only  the  zeroth  and  the  first- 

order  terms  of  the  series,  it  may  be  verified  that  the  zeroth  order  term  is 

equal  to  the  square  of  the  expected  value  E[c  ]  of  r  .  The  first-order 

b  b 

approximate  autocovariance  function 

K^(t)  =  Rfe ( r )  -  E2  T  ]  (27) 

is  therefore  a  linear  function  of  r.-Ci)  ,  rf~^(t)  ,  and  rf.^(i)  the 

12  12  2  12 

Fourier  transforms  of  which  are,  respectively,  S  ( to)  ,  -w  S(w)  ,  and 

4 

w  S(oj)  .  Thus,  by  taking  the  Fourier  transform  of  Equation  27,  we  have  the 
approximate  spectrum  of  the  breaking  wave  simply  related  to  S(ui)  as 


in  which 


Sb(u)  =  F(w)  S (u) 


(28) 


(29) 


is  a  fourth  order  polvnominal  function  of  w  and  may  be  looked  upon  as  a 
filter  function  which  accounts  for  the  effects  of  wave  breaking  on  the 
spectrum  S(w)  of  the  ideal  waves. 

21.  In  Equation  29, 


A1 

(4) 

r 

A„ 

(2) 

2 

r 

(  30) 


Aj  =  BN  +  Q(-B)  >  0 


(31) 


12 


and 


22,  To  examine  the  properties  of  the  filter  function  F(uj)  and  the 

breaking  wave  spectrum  SL (w)  ,  it  is  first  noted  that  the  value  of  6  ,  which 

b 


is  the  ratio  of  Kg  and  the  standard  deviation 


of  the  surface  accel¬ 


eration  of  the  ideal  waves,  may  be  given  a  rough  estimate.  By  referring  to 
Figure  1,  let  us  assume  that  the  acceleration  in  those  portions  of  the  sur¬ 
face,  where  j £ |  reaches  or  exceeds  Kg  ,  remains  at  the  value  of  Kg  ;  but 
in  the  remaining  portion  of  the  surface  the  acceleration  vanishes.  The  stan¬ 


dard  deviation 


is  therefore  equal  to  KgfA  /A  ,  and 

D 


f  =  1 / *  A  /A  where  A„  and  A  are,  respectively,  the  area  of  wave  surface 

n  d 

with  r,‘  >  Kg  and  the  total  area.  The  ratio  A  /A  is  normally  a  small 

D 

quantity  so  that  8  may  he  expected  to  be  larger  than  unity.  For  example,  in 
stornv  situations,  the  ratio  A  /A  mav  be  as  high  as  1/4  giving  8  =  2  ; 

n 

whereas  in  calmer  situations,  if  the  ratio  A  /A  is  equal  to  1/9,  r  is 

n 

approximately  equal  to  3. 

23.  Having  established  that  ~  is  larger  than  unity,  by  employing 
variously  the  series  representation  and  thp  asymptotic  behavior  of  Qf*')  for 
large  values  of  its  argument  (Abramowitz  and  Stegun  19681,  it  may  be  verified 


that  A.,  >  0  and  A  /A.,  >>  i 


S  ince 


(4)/(21 
r  /  v 


the  latter  quantity  is  in  fact  the  characteristic  wave  frequency  w  (see 

Equation  3),  it  is  seen  that  >>  ^  in  view  of  Equation  30.  The 

filter  function  F(v)  is  a  monotonical 1 v  decreasing  function  'f  w  for 

2 

0  ■  ■■  <"  Wj  decreasing  from  F(o)  =  A^  to  Ffw^)  =  0  .  Beyond  us  =  , 

F(wl  increases  indefinitely.  The  range  of  frequency  of  wind  waves  of 
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practical  interest,  however,  is  usually  limited  to  within  0  <  uj  <  (>>  w) 

as  the  numerical  results  to  be  presented  later  will  show.  The  manner  in  which 

F(w)  varies  with  6  ,  a  measure  of  the  sea  state,  may  be  seen  by  taking  the 

derivative  of  with  respect  to  8  .  It  may  be  verified  that  is  a 

monotonically  increasing  function  of  8  ,  and  A^  approaches  unity  as  8 

approaches  infinity  which  means  that  in  mild  seas  6  and  are  both  rather 

large  so  that  F(oj)  =  1  for  0  <  w  <  to,  and  S,  (oj)  =  S(w)  .  No  wave 

]  b 

breaking  takes  place,  and  the  original  ideal  wave  spectrum  remains  unchanged. 
In  high  seas,  on  the  other  hand,  Aj  <  1  and  so  is  F(oj)  for  0  <  u  <  w  ^  . 
Thus,  the  original  wave  spectrum  is  reduced  as  a  consequence  of  wave  breaking, 
as  expected. 


PART  IV:  WAVE-CURRENT  INTERACTIONS 


24.  Consider  a  unidirectional  linear  wave  train  entering  a  region  of 
current.  Let  the  current  be  steady  in  time  and  the  flow  velocity  U  ,  con¬ 
sidered  positive  in  the  direction  of  the  waves,  be  uniformly  distributed  in 
the  vertical  direction.  For  each  wave  component,  the  apparent  frequency  w 
in  a  stationary  frame  of  reference,  is  related  to  the  relative  or  intrinsic 
frequency  in  the  frame  of  reference  moving  with  the  current  as 


u)  —  w  +  k  U 

a  r 


(34) 


where  *  and  the  wave  number  k  are  related  as 
r 


2 

w“  =  gk  tanh  kd 
r 


(35) 


ing  the 

energy 

flux  bal 

ance 

of  wave 

act  ion 

(Hedges , 

Bur 

spectrum 

S(cu  ) 

a 

,  under 

the 

spectrum 

in  qui 

escent  d 

eep  ■ 

and 


SC-  )  =  77-r^2 - 11  S  (w  ) 

a  l  +  c  w  o  a 
gr  a 


(36) 


wh  ere 


go 


2  „ 

a 


(37) 


and 


gr 


4  1  + 


2kd 


sinh  2kd  /  k 


(38) 


The  subscript  "o"  is  used  to  refer  to  quantities  evaluat 
zero  current  condition. 

26.  In  the  relative  frame  of  reference,  the  wave 
be  obtained  from  S(w  )  in  Equation  36  by  changing  the 
(Hedges,  Burrows,  and  Mason  1979)  as 


ed  in  deep  wa 

spectrum  S  ( - 
frame  of  refe 


ter  in 

)  mav 
r 

rence 
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(39) 


dijJ 

S («»  )  =  S(w  )  — ^ 
r  a  du> 


27.  The  above  Is  but  a  brief  exposition  of  the  basic  equations  for  the 

determination  of  the  ideal  wave  spectra  S(w  )  and  S(m  )  for  waves  in  water 

a  r 

of  finite  depth  in  the  presence  of  current.  Details  of  many  of  the  considera¬ 
tions  and  operations  involved  are  well  explained  in  Hedges,  Burrows,  and  Mason 
(1979).  For  example,  the  report  gives  an  account  of  the  solutions  of  Equa¬ 


tions  39  and  3 S  given  the  values  of 


and  r 


frequency  of  (and 


^  )  in  a  negative  current  (c 
a  gr 


discusses  the  cutoff 
!  U  ! ) ,  and  shows  a 

numerical  scheme  by  which  the  transformation  of  the  spectrum  from  the  station¬ 
ary  to  the  relative  frame  of  reference  and  vice  versa  may  be  achieved. 

28.  To  account  for  wave  breaking,  the  spectrum  S(w  )  given  by  Equa¬ 
tion  39  mav  be  used  as  the  original  ideal  wave  spectrum  in  place  of  S  ( wl 

1) 

.2 


(2)  (41 

in  Equations  I S ,  16,  and  17  t or  the  calculation  of  r,  r  and  r  from 


which  the  mean  value  E[4^)  tmd  the  niean-square  value  E  |^J  °f  the  breaking 
waves  are  obtained  from  Equations  14  and  24,  respectively.  Similarly, 

the  breaking  wave  spectrum  in  the  relative  frame  of  reference,  denoted  by 

)  may  be  obtained  from  Equation  28  with  the  ideal  wave  spectrum  S(w) 

replaced  by  S(w^)  in  Equation  39.  Finally,  the  breaking  wave  spectrum  in 
the  stationary  frame  of  reference  is  determined  from 


S,(v  )  =  S  U  J  ~ 
b  a  hr  dw 


(40) 


by  changing  the  frame  of  reference. 
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PART  V:  NUMERICAL  RESULTS 


29.  The  preceding  development  enables  a  consideration  of  the  effect  of 
wave  breaking  on  the  mean  value,  the  mean-square  value,  and  the  spectrum  of  a 
unidirectional  deepwater  wave  train  that  is  free  of  current,  propagating  over 
a  gently  varying  sea  bottom  with  straight  and  parallel  contours  normally 
incident  toward  a  straight  shoreline,  where  it  meets  an  adverse  horizontal 
variable  current  steady  in  time  and  uniformly  distributed  with  depth.  The 
following  computation,  though  not  entirely  realistic,  treats  the  current  speed 
as  a  constant. 

30.  Let  the  deepwater  wave  spectrum  be  the  Wallops  spectrum 
(Huang  et  al.  1981)  which  takes  the  form 


S  (u) 
o 


2 

qg 

m  5-m 

0)  GJ 


exp 


(41) 


where 

a  =  coefficient  defined  in  Equation  44 
u!o  =  parameter  of  Wallops  wave  spectrum 
The  quantity  m  gives  the  magnitude  of  the  slope  of  the  spectrum  (on  log-log 
scale)  in  high  frequency  range  and  is  given  by 


m 


log(2T,V) 

log  2 


(42) 


where 


§ 


(43) 


is  the  significant  slope  of  the  waves,  >  =  2r/k  being  the  characteristic 

o  o 

wave  length.  The  quantity  f  is  given  by 
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m-1  „ 

m  A _ (2Ti§r 

,  m-5  Pf(m-1)‘ 

“  rL~ 4“ . 


(AA) 


where  T(*)  is  the  Gamma  function  (Abramowitz  and  Stegun  1968).  The  Wallops 
spectrum,  therefore,  is  seen  to  depend  on  two  parameters,  §  and  w  ,  the 
frequency  corresponding  to  the  peak  of  the  "single-peak"  Wallops  spectrum. 

31.  For  current  speed  of  U  =  -2m/sec  ,  §  =  0.015  (the  value  of  § 

rarelv  exceeding  0.025  in  the  field),  and  uj  =  0.6  rad/sec  ,  the  quantities 

r  2  n  ° 

E[f;,  ]  ,  E  ,  S.  (to  )  ,  and  S,  (w  )  are  computed  for  various  values  of 

b  L  °J  b  r  b  a 

water  depth  d  .  The  solutions  are  carried  out  in  an  iterative  manner;  that 

is,  upon  obtaining  S  (u  )  in  Equation  28,  it  is  treated  as  the  original 

b  r 

ideal  wave  spectrum,  and  the  solution  process  is  repeated  until  convergence  is 

reached.  Based  on  the  final  values  of  S,  (w  )  ,  the  quantities  E  [  r  ]  , 

2  hr  h 

E  t,  ,  and  S.  (w  )  are  then  determined.  The  results  presented  in  the 
b  b  a 

following  are  obtained  after  four  cycles  of  iteration.  It  should  be  mentioned 

that  the  above  iterative  scheme,  strictly  speaking,  is  not  valid  since  some  of 

the  assumptions  underlying  the  derivation  of  these  quantities  are  violated 

because  it  was  original Iv  assumed  that  the  ideal  waves  must  be  zero  mean  and 

Gaussian.  Our  results  show,  however,  that  the  mean  value  of  c,  is  insig- 

b 

nificnntly  small,  and  preliminary  investigation  indicates  that  the  breaking 

wave  elevation  c,  deviates  but  slightlv  from  Gaussian, 
b 

32.  In  Figure  2,  Ef^l  is  plotted  as  a  function  of  k  d  for  k^d 

ranging  between.  3  ard  0.5  where  k^  =  i.*/g  is  the  characteristic  deepwater 

wave  number,  c  being  given  bv  Equation  3  with  S(.  )  replaced  by  S  (  ..)  , 

-  -  i/o  _  ° 

the  Wallops  spectrum.  If  we  denote  by  k  =  kytanh  k^d  according  to 

Fquation  8,  these  values  of  k  d  correspond  to  kd  1  A.l  and  1  (or  d  = 

o 

81.3  and  13.6  ).  It  is  seen  that  F  f  -  ^ 1  is  always  negative,  as  expected, 
and  indeed  very  small. 

33.  Figure  2  also  gives  the  standard  deviation  >r  (see  Equation  151 


the  elevation  of  the  ideal  waves  and  ►'r. 


rF; 

r 

9 

-  f  rt  ] 

v  ' 

[’bj 

1  V 

breaking  waves.  While 


decreases  monotonical lv  from  deep 


that  of  the 
to  shallow 


water,  >r  first  decreases  slightlv.  Bevond  k  d  =  0.6  shoreward,  however, 

o 

it  begins  to  rise  because  of  shoaling.  Owing  to  the  relativelv  small  value  of 
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Figure  2.  Mean  value  E [ ^  ]  of  breaking  waves,  standard  deviations  fr 

of  ideal  waves,  and  t'r  of  breaking  waves  plotted  against  k  d  for 

b  ° 

§  =  0.015  ,  -JJ  =  0.6  rad/sec,  and  0  =  -2  m/sec 
o 

5  and  the  strong  negative  current  speed  used,  wave  breaking  is  seen  to  occur 
everywhere  but  more  so  in  shallower  than  in  deeper  water. 

34.  For  k  d  =  3,  2,  1  and  0.5  ,  the  spectra  S  (v  )  ,  S(w  )  ,  and 

o  r  r  r 

g  ( .  )  are  in  Figures  1,  5,  7,  and  9,  respectively,  and  those  in  the 
h  '  r 

stationary  frame  of  reference  S  ( w  )  ,  S(u  )  ,  and  S.  (a  )  are  given  in 

o  a  a  n  a 

Figures  4,  6,  8,  and  10,  respectively.  An  adverse  current  feeds  energy  into 
the  wave  system  so  that  the  ideal  wave  spectra  always  exceed  those  in  deep 
water  where  there  is  no  current.  Wave  breaking,  however,  dissipates  wave 
energy,  and  the  breaking  wave  spectra  are  seen  to  fall  below  S  (O  •  ('lose 
examination  of  these  spectra  also  shows  that  this  pattern  ot  variation  with 
water  depth  is  consistent  with  that  of  the  standard  dev  :  -a  t  ions  shown  in 
Figure  2 . 

35.  As  mentioned  earlier,  the  quantity  c  g  i  v*  s  an  in.!  i  eat  ion  .-•!  the 
extent  of  wave  breaking  and  is  expected  to  he  larger  than  urit  \  ,  with  t  he 
larger  values  corresponding  to  milder  sea  state.  It  was  also  shown  that  t be 
quantity  w  ,  the  cutoff  frequency  of  the  breaking  wave  spectrum  Mw  i  r. 
Equation  30,  is  expected  to  be  much  larger  than  ,,  ,  the  character  1st  ii  wave 
frequency.  It  is,  therefore,  of  interest  to  exnnim  the  variation  d  t  hr*-.- 
two  quantities  as  the  waves  move  toward  the  sliore.  in  Figure  l.1  the 


'  r  ad  .'sec! 


Figure  3.  Deepwater  wave  spectrum  ,  idea 

wave  spectrum  S/ou^)  ,  and  breaking  wave  spectrum 

S,  (w  )  in  relative  frame  of  reference  for  k  d 
b  r  o 

=  3,0  ,  §  =  0.015  ,  oj  =  0.6  rad/sec,  and 

o 

U  =  -2  m/sec 


0  «  OS  3  •  3  ^  ;l  39 


Figure  4.  Deepwater  wave  spectrum 

S  (ui  )  ideal  wave  spectrum  S(w  )  , 
o  a  a 

and  breaking  wave  spectrum  ) 

in  stationary  frame  of  reference  for 

k  d  =  3.0  ,  §  =  0.015  ,  u 

o  o 

=  0.6  rad /sec,  and  U  =  -2  m/s ec 


Vi  I 


Figure  5.  Deepwater  wave  spectrum  S  (~  )  ,  idea 

wa.e  spectrum  S(i^)  ,  and  breaking  wave  spectrum 

.  t.  ~  )  in  relative  frame  of  reference  for  k  d 
hr  o 

-  2.0  ,  §  =  0.015  ,  w  =0.6  rad /sec,  and 

o 

!'  =  -2  m/sec 


Figure  6.  Deepwater  wave  spectrum. 

S  ( *.  )  ,  ideal  wave  spectrum  S(  .. 

o  a  a 

and  breaking  wave  spectrum  1 

in  stationary  frame  of  reference  for 


Ki  mire  .•  .  Deepwater  wave  spectrum  Si  «;  )  .  tciea 

wave  spectrum  Sfw  )  ,  and  breaking  wave  spectrum 

S,  <  .  )  in  relative  frame  of  reference  for  k  d 

h  r  o 

=  l.(i  ,  <>  =  0.015  ,  vj  =0.6  rad /sec,  and 

o 

1'  =  -2  m/ sec 


Figure  8.  Deepwater  wave  spectrum 

.si.  )  ,  ideal  wave  spectrum  S(  ,  )  , 

n  a  a 

and  breaking  wave  spectrum  S^(  > 

in  stntinnarv  frame  ot  reference  lor 
kd=  1.0.  §=0.015,  c 


•J.  I 
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Figure  11.  Wave  breaking  parameter 
wave  spectrum  plotted  against  k^d 

V  = 


8  and  cutoff  frequency 

for  §  =  0.015  ,  u<  =  0.6 
o 

-2  m/sec 


oi.  of  breaking 
rad /sec  ,  and 


quantities  8  and  are  plotted  against  k^d  .  The  value  of  6  decreases 

from  1.8  in  deep  water  to  about  1.4  at  k  d  =  0.5  and  that  of  w,  from 

o  1 

b.O  to  4.0  rad/sec. 

36.  Although  our  primary  objective  in  this  study  is  to  devise  a  method 

for  the  calculation  of  the  breaking  wave  spectrum  under  the  influence  of  an 

adverse  current,  the  method,  as  is  obvious,  may  be  applied  to  the  special  case 

in  w'hich  there  is  no  current.  The  results  presented  in  the  following  are  for 

1’  =  0  ,  §  =  0.015  ,  and  w  =  0.6  rad/sec  for  k  d  =  3  to  0.5  correspond i ng 

o  o 

to  kd  =  3  to  0.74  . 

37.  In  Figure  12,  F  [  t; ,  ]  ,  *' r  ,  and  *  r,  are  given  as  functions  of 

h  b 

k^d  .  As  expected,  Eft^)  is  always  negative  and  oven  smaller  than  when 
F  =  -2  m/sec  in  Figure  2.  Because  of  the  small  value  of  significant  slope 


5  =  0.015  used,  the  two  curves  >r  and  rr^  are  practically  indistin¬ 


guishable  until  k  cl 
o 


-  1  .5  (i 


kd  <  1.6,  d  <  41.7  j  when  wave  breaking  becomes 
noticeable  and  they  begin  to  diverge  from  each  other. 

38.  The  variation  of  the  quantities  8  and 
Figure  13.  The  value  of  8  changes  from  2.9  to  1.9,  and  that  of 


with  k  d  is  shown  in 
o 


1 


from 


20  rad/sec  to  3  rad/sec  as  k  d  goes  from  3  to  0.5.  Comparison  of  Figure  13 


with  Figure  11  for  the  case  of  l!  =  -2  m/sec  shows  that  the  quantities  f 
and  uu  undergo  more  variation  in  the  present  case  of  1=0  that:  in  the  case 
l'  =  -2  m/sec  because  in  the  latter  case  the  current  dominates  the  flow  field 
so  that  wave  breaking  is  uniformly  present  regardless  of  the  locale  under 
considerat ion . 

39.  The  spectra  of  the  deepwater  waves  S  (w)  ,  the  local  ideal  waves 

S(oj)  ,  and  the  breaking  waves  S^('w)  will  not  he  shown;  but  the  peak  values 
(in  square  metres  per  second)  are  recorded  in  the  following  tabulation: 


PART  V  I  : 


CONCLUSION 


40.  We  have  given  in  this  report  a  method  to  compute  wave  spectra  in 
water  of  finite  depth  taking  into  account  the  effect  of  wave  breaking  and  con 
sidering  the  presence  of  current.  The  breaking  wave  spectrum  is  simply 
related  to  that  of  the  original  ideal  waves.  The  method,  however,  is  approxi 
mate  because  the  following  is  assumed:  (a)  there  exists  an  ideal  original 
wave  train  which  is  linear  and  Gaussian;  (b)  the  wave  breaking  model  is 
heuristic,  and  some  approximations  have  been  introduced;  and  (c)  the  higher 
order  terms  in  the  expression  of  the  breaking  wave  spectrum  are  ignored.  As 
such,  the  model  should  only  he  applied  to  the  energy  containing  part  of  the 
spectrum  and  must  not  he  used  for  points  too  close  to  the  shore  where  all  the 
assumptions  underlying  this  model  will  be  violated. 
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APPEND  I X  A:  DERIVATION  OF  E  f  r. 


1.  To  obtain  E ( r j ' nH''+H"+ ]  ,  the  concept  of  conditional  probability 
and  conditional  expectation  to  reduce  the  number  of  random  variables  is  used, 
thus  (Papoulis  1965),* 


F^"l  ’2Hl'+H2+I 


•  '2]1 


( A1  ) 


where 


(Ad) 


is  the  conditional  expected  value  of 


g  i  v  e  n  '  a  n  d 


and 


O  1 


(,  A 1 1 


is  the  jointly  Gaussian  conditional  probability  density  function  of  '  and 

'  ,  given  .  and  r  where 
1  1 

t  =  conditional  covariance  coetficient  function  of  r  ^  and  r,„ 

»  O 

■  ’  j  ,  ' ,,  =  conditional  variance  functions  of  ’  ^  and  r. 

M|  ,  p =  conditional  mean  value  functions  of  '  and 

The  five  parameters,  o  j  0  ,  a~  ,  ,  p^  ,  and  u ,  ,  may  all  be 

determined  using  the  linear  mean-square  estimation  technique  (Papoulis  1965); 
that  is. 


*  References  cited  in  the  Appendix  can  he  found  in  the  References  at  the  end 
of  the  main  text. 


A! 


3.  In  the  above,  r  =  r  ,(o)  ,  r  =  ,  (nj  ,  anil  r  ■  t ,  ,  -i'1 

are  Riven  in  terms  of  S  (.  i  ns  indicated  in  Iquutions  1  lb,  and  I  ,  reaper 

/  (  2  i  ,  (  •.  i  ,  \ 

tively.  The  argument  t  l  the  time  lag  in  r  ^  f  •  '  ,  r^,  i)  ,  and  r  ,  ,  1  ’  1  ) 

is  omitted  for  brevity.  The  quantities  Uj  ■  u.,  .  '|.  ,  •  -md  i  ^  ,  an 


A 


( 2 )  (4) 

all  functions  ot  r  0(t)  .  rj^  (t)  ,  and  rj?  (t)  and  hence  are  functions  of 

4.  The  conditional  expected  value  E  f  4  |  i;  ,  £ ^  1  is  seen  to  be  the 

conditional  correlation  function  of  £  and  and  is  therefore  by  defini¬ 

tion  given  bv 

K [ r , | C 1 » to  1  =  U  j  U  2  +  Pi2°l°? 


/  2  2  ..  ■■  2  .2  \ 

'  al  +  bl)ClC2  +  albl(Cl  +  ^2  )  +  p  1 2°  1  2 


The  expected  value  sought  is,  therefore,  from  liquation  Al  , 


=  ff  [(ai +  bi }r- 1  ^2 +  aibi(’i +  '&) +  °  1 2 °  1  ° 2 


f  ••  v  (i;,  ♦  C ) d t  d r' , 


whore  :  ■  ( 4 .  ,  r )  is  the  jointly  Gaussian  probability  density  function 

’  1  ’  r  2 

1!  the  zero-mean  random  variables  r  and  ;jrj  whose  variances  are 

r .  1  r.pi  -  ( 4  j 

K;.  :  =  Mi  =  r  ,  and  whose  correlation  coefficient  function  is  p  '(;) 
L.  '  J  L  “  J  ‘  *- 


<  4i  (41  /  (41 

r  =  r  .  ,  /  r 


The  above  integrals  rav  all  be  carried  out  giving 


2H'l'+H2+  I 


=  +  SjFl  +  2aihlF2  4  rU  VjF1 


v  ^  '  1 

r  ,  =  r  <1 


.1  1  /  " 


141  '  (4) 

Z(  e)Z(;:)  -  2pc  ,  7.!  pH).  )  +  c '  F.,  V  Al  Al 
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(4) 
1  2 


Z(B)Z(fi) 


( 4 )  V 

0  1  2  j 


Z(£)0(S-)  +  F. 


(A1  5) 


and 


F,  -  l-(-C 

(pee  Abramowitz  and  Stegun  1968).  Here, 
omi t  ted 


-t  , 


the  argument 


1  - 


I  + 


(4) 
"12 
.  (4) 

‘  1  2 


1/2 


(AIM 


(  4  )  ,  . 

in  P . 0  ( t '  is 


(  A1  .  i 


and  r  ,  '/.(")  ,  and  Of*''  are  defined  in  Equations  19,  20,  and  22, 

respectively. 

6.  The  expected  value  in  Equation  A!  is  a  nonlinear  function  of 
r  j i  ,  Tjt,  (:)  ,  and  r  ^ (  T  )  and  rrav  be  expanded  by  the  Taylor  series. 
Bv  retaining  on  1 v  the  zeroth  and  the  first  order  terms,  it  is  given 
approx ima t el y  by 


r  r 
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:  =  a,  >- '  *  '  7."  (t)  +  a  T  r  (  ^  c  !  Z,  [  - -Z  (  r  )  +  0(-B)  1 ' 


( A 1  8  i 


(4) 


+  2a  hr  Q(-t)  [-(-/if:: )  +  f  ,  a  '.,0"  ( -r ' 


wh  c  r  o 


( 4 ) 


[  A  1  <»  1 


A  4 


(2>  (.+  )  (i)  (?) 

_rl ?  r  rI?  J  r 


(2)  ...  (2) 

„  (?)  r _  ( 4 )  r 

ri?  (4)  r 1 2  (4) 

r  r 


The  integration  in  Fquatiori  A 1 2  may  also  he  facilitated  by  employing 
the  Hermit*?  pel vroirin.il  series  representation  CKrdelv  et  al.  1953)  as  follows: 


(x~  +  v“  -  2pxvl 


2 1 1  -  p" ) 


=  2~  ^  ^  pnhn  (x)  (y )  7.  f  x)  Z  (y) 


h  :<)  = 


,n„  ■> 

d  /.  ( y.  1 

nn  ,  n 
,  -  1  )  dx 

,  J  /  ?  ThO 


is  the  liernite  polvnominal  function.  Upon  expanding  the  sointlv  i-aussian 

pr  >hah  i  1  itv  density  function  fy  •;  into  the  dermite  ‘•cries 

’  1  ’ 

represent, it  ion,  if  is  seen  that  the  integrals  in  Fquat  ion  A 1 ?  may  he  carried 
our  easily.  Kv  retaining  only  the  terms  involving  n  =  0  and  1  in  the 
series,  Kquation  A18  mav  be  obtained. 
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APPENiJ  1 X  ll: 

Total  wave  surface  area  and  wave  surf are  area  with 
"h!  •  Kg 

Puant  ities  del  ined  in  Kquations  31  and  32, 
reaper t i ve 1 v 

.Amplitudes  of  ideal  and  breaking  waves,  respective  1  \ 

Quantities  del  ined  in  Equations  A6  and  .47  whose 
approximate  values  are  given  in  Equations  AIQ  and  A70 

Quantities  defined  in  Equations  A7  and  A6  whose 
approximate  values  are  giver  in  Equations  A20  and  A!° 

Oruip  velocities  defined  in  Equations  37  and  38, 
respec t  i  ve 1 v 

Parameters  used  in  Equation  21 
Water  depth 

Expected  value  of  the  quantity  enclosed  in  the 
brackets 

Conditional  expected  value 

Exponential  integral  defined  in  Equation  23 

Quantities  defined  in  Equations  .414,  A15,  and  Alb 

Filter  function  defined  in  Equation  29 

Joint  probability  density  function  or  the  Gaussian 
zero-mean  random  variables  •’  and 

Conditional  point  probahilitv  density  function  of 
r.  j  and  '  ,  given  ^  and 

Cravj  tat  i  ona  1  acce  lerat  i  or. 

Heaviside  unit  step  function 

Abbreviations  for  H(’l  ,  \U  r  4  Kg!  ,  and 
H(-r,  -  Kg)  ,  respectively 

llermite  polynomial  function  defined  in  Equation  A23 

Coefficient  defined  in  Foliation  I  5 

Aut  ('covariance  function  oi  breaking  wave  elevation 


Wave  number  and  characteristic  wave  number,  respec¬ 
tively,  in  water  of  finite  depth 

Wave  number  and  cliarac  t  er  1  st  i  c  wave  number,  respec¬ 
tively,  in  deep  water 

Probe  ili tv  function  defined  in  Equation  21 
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r  ,  r 
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r  , 
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■.it', 
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ii 


M  *  gn  i  t  mil'  o'  tilt'  ■  !  ope  of  the  i-.’.'i  1  !  ■  ! 
doepwat  or  vavi'S  in  the  high  freouercv  r. 
i  in;  scale  given  in  liquation  4 J 

Quantities  del  ineJ  in  !'<|tiat  ions  .'4  usd 
respect i ve ! v 

I'ummv  index 

Parameters  used  in  Kqunt i nns  ?1  and  A.  1’ 
I’rohahilirv  function  defined  in  Fount  im 
An  t  t'cor  re  1  a  t  i  on  function  of  breaking  v:,v 

V'1 

Quantities  defined  in  Kquntions  !a,  lb, 
Variance  of  breaking  wave  elevation  " 


Ci’  rr  e  i  a  t  i  or,  functions  K  [  ■  ■  ■  , 


r  i 


!■  j  '  .  ’  ,  I  of  the  ideal  waves,  respect  had 
Ideal  and  breaking  wave  spectra,  respect 
V.'ave  spectrin-  in  deep  water 


Ideal  and  breaking  wave  spectra,  respect 
re  i.at  i vt  *  rame  of  reference  it;  the  prose 
curron t 

Tim 

T  i  tto  instants  t  4  -  and  t  ,  -osrua  t  i\a 
Current  speed 

Quantity  defined  in  fount  ion  hi 
Qunrr.v  variables 

Probabilitv  function,  defined  in  1  quarter 

Coefficient  defined  in  liquation 

Wave  breaking  parameter  defined  in  Fount 

The  Carina  funct  ion 

Ouantitv  defined  in  equation  AM 

'•'pretral  bandwidth  parameter  defined  in  1 

Illevat  inns  id  ideal  and  breaking  wave  , 

Characteristic  wave  length 

Conditional  mean  value  functions  of  • 
respec  t  i  vel  y  ,  given  and  c 

Conditional  covariance  coefficient  funct 

and  r ,,  given  '  ^  and  ’  in  Kquatinn 
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'.aw t 


i  rn  « '  t 
Ain 


H2 
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u) 

CO 


(U 
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Subscripts 

1,  2 


o 


Correlation  coefficient  function  of  r  and 

Conditional  standard  deviation  of  4  and  , 

respectively,  given  £  and  c. 

Time  lag 

Quantity  defined  in  Equation  A17 
Wave  frequency 

Characteristic  wave  frequency  defined  in  Equation  3 

Wave  frequency  in  stationary  and  relative  frames  of 
reference,  respectively 

Parameter  of  Wallops  wave  spectrum 

Cutoff  frequency  of  breaking  waves  given  in 
Equation  30 


Quantities  evaluated  at  time  instants  tj  and  t0  , 
respectively 

Quantities  evaluated  in  deep  water  in  7ero  current 
condition 


Svmbols 


§ 


Differentiation  with  respect  to  time 
Significant  wave  slope  defined  in  Equation  43 


B3 


